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GAP PROBABILITIES FOR THE CARDINAL SINE 

JORGE ANTEZANA, JEREMIAH BUCKLEY, JORDI MARZO 
AND JAN-FREDRIK OLSEN 

Abstract. We study the zero set of random analytic functions generated by a 
sum of the cardinal sine functions that form an orthogonal basis for the Paley- 
Wiener space. As a model case, we consider real-valued Gaussian coefficients. 
It is shown that the asymptotic probability that there is no zero in a bounded 
interval decays exponentially as a function of the length. 



A simple point process in R is a random integer-valued positive Radon measure 
on R that almost surely assigns at most measure 1 to singletons. Simple point 
processes can be identified with random discrete subsets of R. In this paper, we 
study 'gap probabilities' of the simple point process in R given by the zeros of the 
random function 



where a n are i.i.d. random variables with zero mean and unit variance. Kol- 
mogorov's inequality shows that this sum is almost surely pointwise convergent. 
In fact, since 



converges uniformly on compact subsets of the plane, this series almost surely 
defines an entire function. If we take a n to be Gaussian random variables then / 
is a Gaussian analytic function (GAF) (see |HKPV[ Lemma 2.2.3] for details). 

We shall be chiefly concerned with the functions given by taking a n to be real 
Gaussian random variables. We denote by rif the counting measure on the set of 
zeros of /. These functions are an example of a stationary symmetric GAF and 
Feldheim jF] has shown that the density of zeros is given by 



where z = x + iy, m denotes the planar Lebesgue measure, \i is the singular 
measure with respect to m supported on R and identical to Lebesgue measure 
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(Here S is defined only for y ^ 0, in fact the atom appearing in Q is the distri- 
butional derivative at 0.) We observe that since S(y) = 0(y) as y approaches zero 
there are almost surely zeros on the real line, but that they are sparse close by. 
Moreover the zero set is on average uniformly distributed on the real line. We are 
interested in the 'gap probability', that is the probability that there are no zeros 
in a large interval on the line. Our result is the following asymptotic estimate. 

Theorem 1. Let f be the symmetric GAF given by the almost surely convergent 
series 

sin tt(z — n) 



ir(z — n) 



where a n are i.i.d. real Gaussian variables with mean and variance 1. Then, there 
exist constants c, C > such that for all r > 1, 

e- cr <P(#(Z(/)n(-r,r)) = 0)<e- Cr . 

Remark 1. If instead of considering intervals we consider the rectangle D r = 
(— r, r) x (—a, a) for some fixed a > 0, then we compute a similar exponential 
decay for P(#(Z(/) n D r ) = 0) 

Remark 2. Suppose that the a n are i.i.d. Rademacher distributed. I.e., each a n 
is equal to either —1 or 1 with equal probability. Since f(n) = a n for n G N, it 
follows that if not all a n for |n| < iV are of equal sign, then by the mean value 
theorem, / has to have a zero in (-N, N). As is shown below, the remaining two 
choices of the a n for \n\ < N each yield an / without zeroes in (-N, N). Clearly, 
this gives a probability of 2(1/2) 27V for / to be without zeroes there. 

Remark 3. The Cauchy distribution is given by the density 

P{x) = ~ 2 ~ r - 

7T X 2 + 1 

Whereas the Rademacher distribution is in some sense a simplified Gaussian, the 
Cauchy distribution is very different: It has neither an expectation, nor a standard 
deviation. If we suppose that the a n are i.i.d. Cauchy distributed, it is not hard 
to see that with probability one the sum J^a n /n diverges, whence the related 
random function diverges everywhere. 

The main motivation for our work comes from the 'hole theorem' proved by 
Sodin and Tsirelson in |ST] for point-processes uniformly distributed in the plane. 
The authors consider the GAF defined by 
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where a n are i.i.d. standard complex normal variables. For this function, the 
density of zeros is proportional to the planar Lebesgue measure and the authors 
compute the asymptotic probability that there are no zeros in a disc of radius r 
to be e~ cr , where c > 0. The analogy with the GAF we consider becomes clear 
if we note that (^=y)^L constitutes an orthonormal basis for the Bargmann-Fock 
space 

T = {/ G H(C) :\\f\\%= I |/(z)|V 2 N 2 ^ < +00}, 

Jc 71 

where m is the planar Lebesgue measure. We are replacing these functions with the 
sine functions, which constitute a basis for the Paley- Wiener space. An important 
caveat is that, though / is constructed from an orthonormal basis, / is almost 
surely not in the Paley- Wiener space, since the sequence of coefficients a n is almost 
surely not in £ 2 (Z). It is not hard to see however that / belongs almost surely to 
the Cartwright class. 



2. Proof of Theorem [TJ 

2.1. Upper bound. We want to compute the probability of an event that con- 
tains the event of not having any zeroes on (-N, N), for iVeN. One such event 
is that the values f(n) have the same sign for |n| < N. The probability of this 
event is 



F(a n > for |n| < A" or a n < for|n| < Nj = 2(1/2) 2N+1 = e~ CN , 
for some constant C > 0. 

Remark 4. The same upper bound holds when a n are i.i.d. random variables with 
< P(a n > 0) < 1 for which the random function Ylnez OnSinc(:r — n) converges. 

2.2. Lower bound. To compute the lower hole probability, we use the following 
scheme. First, we introduce the deterministic function 

2N 

fo{ x ) — /J sinc(x — n). 

n=-2N 

We show in Lemma [T] that it has no zeroes on (—N,N), and we find an explicit 
lower bound on (-N, N) for it. This lower bound does not depend on N. Second, 
we consider the functions 

2N 

fi(x)= (a n — 1) sine (a; — n) and f2(x)= a n sinc(x — n), 

n=-2N \n\>2N 

which induce the splitting 

/ = fo + h + h- 
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We show that for all x G [— N, N] we have |/i(x)| < e with probability at least 
e -cN f or } ar g e ]\f anc i some constant c > 0. Moreover, we show that 

P( sup |/ 2 (ar)|<e 

\xe[-N,N] 

is larger than, say, 1/2 for big enough N. As the events on f\ and f 2 are clearly 
independent, the lower bound now follows by choosing e small enough. 
We turn to the first part of the proof. 

Lemma 1. Given N G N and 

2N 2N , x , n 

fo(x) — 7 sincfa; — n) = sin nx > — r. (2) 

^— ' n{x — n) 

n=-2N n=-2N y ' 

Then, there exists a constant C > such that, for N big enough, 
1 - § < , inf /„(*) < sup f (x) < 1 + 

iV |z|<iV |a;|<JV ^ 

Proof. Let i? = R(N) be the rectangle of length 4iV + 1 and height AN, centered 
at x = 1/2. By the residue theorem, it holds that 

If dz 1 A (-1V 1 



1 E 



2ni Tp (z — x ) sin ttz it ^— ' n — x sin 7rx 
Observe that if we shift around the terms, this yields 

(— l) n sin7rx f dz 

smrcx 



' x — n 2ixi 

n=-2N+l 



R 



(x — sm7r2 

Now, given — iV < x < N, it is easy to bound this last integral by C/N. □ 



Remark 5. The same bound holds for all points 2 in a strip with fixed height 
[— N, N] x [— C, C] for some C > 0. We observe though, that the function fo(z) in 
§2§ is close to zero around Imz = logiV. Indeed, it is smaller than e~ cN there, for 
some c > independent of N. 

2.2.1. The middle terms. Let e > be given, and consider N G N to be fixed. We 
look at the function 

2N sin7rx 2N (-1)™ 

fii x ) = y^A a n~ l)sinc(x -n) = V" (a n - 1) . 

z — ' 7r z — ' x — n 

n=-2N n=-2N 

To simplify the expression, we set b n = (a n — 1)(— l) n . We want to compute a 
lower bound for the probability that, for x G [-N, N], 

\fiW\<e. 

Since, for \n\ < N, we have f\(n) = a n — 1, the condition 

\a n — 1| < e for \n\ < iV 
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is necessary. 

Define B n = 6_2jv + ■ ■ ■ b n for \n\ < 2N with _B_2tv-i 
x ^ Z. With this, summation by parts yields 



27V 

E 



x — n 



27V 

E 



a, 



+ 



B 



-27V -27V 

We now claim that under the event 



x — n)(x — n — 1) x — 2N — 1 



0, and suppose that 
(3) 



27V 



E 



bJ < e, \BJ < e for \n\ < 2N 



(4) 



we have \ f{x) — fo{x)\ < e for \x\ < N, with a bound independent of N. Indeed, 
the second summand at the right hand side of ^ converges almost surely to zero, 
because 

B2N 



X-2N-1 



< 



N' 



Suppose that x e (k, k + 1) and split the first sum in (|3| as 

L>„ \ ^ L) 

T) 4^ 



27V 

E 



(x — — n — 1) ^— ' (x — n)(x — n — l) 

n=-2N V 7V 7 n=fc-l v yv 7 

n^fc-l,fc,fe+l 



Then 



27V 

E 

n=-27V 
n^fc-l,fc,fc+l 



(x — n)(x — n — 1) 



< 



E 



(fc + 1 -n) 2 ' ^ (k-n) 2 ~ 

n>fc+2 v 7 n<fc-2 v 7 



E 



For the remaining terms, the function sinyTX comes into play. E.g., suppose that 
\x — k\ < 1/2, then 



sm7ra 



(x — k)(x — k — 1) 



< 



\x — k — 11 



sin7r(x — A;) 



7r(x — A;) 



< 



e. 



The remaining terms are treated in exactly the same way. 

What remains is to compute the probability of the event E defined by Q. We 
recall that the b n were all defined in terms of the real and independent Gaussian 
variables a n . So the event E above defines a set 



V 



(t- 



27V, 



t 2 N) e 



047V+1 . 



\t n \ <e, |^t n | < e, |n| < 2N 



-27V 



in terms of the values of the a n . Hence, 



F(E) = c 



27V 



< t2 -2N + -4N)/ 2 dt_ 2N ■ ■ ■ dh,K- 



(-27V- 
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. Vn 




FIGURE 1. Illustration of the solid Vr 



AT- 



Here, c is the normalising constant of the one dimensional Gaussian. Since \a n — 
1| < e, it follows that 

F(E) > c 2N e- N ^ 2 J ■■■ J dt. 2N ■ ■ ■ dt 2N = c 2N e- N ^ 2 Vol(V). 

We now seek a lower bound for this euclidean (4iV + l)-volume. 

To simplify notation, we pose this problem as follows. For real variables X\, . . . ,x^, 
we wish to compute the euclidean volume of the solid Vjv defined by \xi\ < e for 
i = 1, . . . , N and 

\x± + x 2 \ < e, 

\xi + x 2 + x 3 \ < e 

\xi + x 2 -\ + x N \ < e. 

One way to do this is as follows. Write y^ = x± + ■ ■ ■ + xn-i, then 

/p / rmm{e,e-y N } \ 

• • - / / dx N J dxi ■ ■ -dxjv-i. 

This is illustrated in Figure [T] Clearly, whenever < 0, the upper limit is e, and 
whenever y^ > 0, the lower limit is — e. Hence, 



Vol(^v) >/•••/ ( / dxpf J dx\ ■ • • dx N ^ 



+ /•■■/ ( / dx N ) dxi ■ ■ • dx^i = eVol(Viv-i) 

J ^Vw_in{?/Ar>0} \J-€ J 



'y ]V -in{?/ J v>o} 
Iterating this, we get 

VoI(Vat) > e*. 
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In conclusion, 

F(E) > e~ cN , 
which concludes this part of the proof. 

2.2.2. The tail. We now turn to the tail term 

fi{ x ) — yj a n sinc(x — n). 

\n\>2N 

Clearly, we need only consider the terms for which n is positive. Set c n = (— l) n a n . 
We factor out the sine factor as above, and apply summation by parts, to get 

L L 



y = y c — 

z — ' X — n z — ' (nr.— 



1 



(x — n)(x — n — 1) 



+ 



Cr 



L- 1 



(5) 



n>2N 2JV+1 

where 

C n = C2N+1 + • • • + C n , C2N = 0. 

We want to take the limit as L — > 00. It is easy to see that the last term almost 
surely tends to zero. Indeed, Cl is a sum of independent normal variables with 
mean and variance 1, and therefore is itself normal with mean and variance 
L — 2N. Moreover, since 



Cr 



x — L — 1 



< 



C, 



N-L 



and the random variable inside of the absolute values on the right-hand side has 
variance L — 2N, it follows that limit is almost surely equal to 0, whence we are 
allowed to let L — > 00 in 

We prove the following. With a positive probability, we have 

L 

cv, 



(x — n)(x — n — 1) 



< e. 



2N+1 

As n 2 ~ \(x — n)(x — n — 1)\ for \x\ < N and n > 2N, it is enough to consider the 
expression 

|Cn| 

n 2 



E 

2JV+1 



The absolute value of a Gaussian random variable has the folded-normal distribu- 
tion. In particular, if X ~ N(0, cr 2 ), then 



E{\X\) = a 



Since, in our case, cr 2 = n — 2N, this yields 



ME 



\Cn\ 



< 



\2N+1 



n- 



E 

2N+1 



Vn-2iV 



< v 

~ ^ (n + 2Nf/ 2 



< 



N 
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Finally, by Chebyshev's inequality, 

/A ICJ \ E(Y) C 

P > Mr < e > 1 ^ > 1 = . 

\ ^ n 2 / e eViV 

\27V+1 / cviv 
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